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ABSTRACT 

We present the GALEX detection of a UV burst at the time of explosion of an optically normal 
Type II-P supernova (PSl-13arp) from the Pan-STARRSl survey at z = 0.1665. The temperature and 
luminosity of the UV burst match the theoretical predictions for shock breakout in a red supergiant, 
but with a duration a factor of ^50 longer than expected. We compare the NUV light curve of 
PSl-13arp to previous GALEX detections of Type IIP SNe, and find clear distinctions that indicate 
that the UV emission is powered by shock breakout, and not by the subsequent cooling envelope 
emission previously detected in these systems. We interpret the ~ Id duration of the UV signal with 
a shock breakout in the wind of a red supergiant with a pre-explosion mass-loss rate of ~ 10~ 3 M 0 
yr^ 1 . This mass-loss rate is enough to prolong the duration of the shock breakout signal, but not 
enough to produce an excess in the optical plateau light curve or narrow emission lines powered by 
circumstellar interaction. This detection of non-standard, potentially episodic high mass-loss in a 
RSG SN progenitor has favorable consequences for the prospects of future wide-held UV surveys to 
detect shock breakout directly in these systems, and provide a sensitive probe of the pre-explosion 
conditions of SN progenitors. 

Subject headings: ultraviolet: general — surveys - supernovae 


1. INTRODUCTION 

The first electromagnetic signature of a core-collapse 
SN is expected to occur when the explosion shock wave 
accelerates through the stellar envelope, and the radi¬ 
ation from the shock escapes and “breaks out” of the 
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stellar surface ( Falk|1978 Klein fc Chevalier|1978 ). This 
radiative precursor to the SN explosion is expected to be 
the most luminous phase of a SN, and yet it is the least 
well observed. This is because the associated burst of 
radiation is brief (1 1000 sec) and hot 10 5 — 10 6 J0) 


( glinnjkove^g I.|2602 Ensman fe Burrows|1992 Klein & 


Chevalier 1978), and precedes the optically-bright phase 


during which SNe are typically discovered. 

These obstacles to discovery were surmounted with a 
bit of luck in a pointed X-ray observation of a galaxy 
in which a Type Ibc SN fortuitously exploded, and a 


berg et al. 

2008 

Modjaz et al. 2009), with a c 

uration and 
ck breakout 
;ar progeni- 

spectral shape that was consistent with sho 
into the thick wind of a Wolf-Rayet (WR) s 

tor ( 

Balberg & Loeb 

2011 

Svirski & Nakar 

2014). The 


presence of a thick wind is expected for a WR star, but 
it was also an important factor in stretching the shock 
breakout timescale long enough (a factor of 100 longer 
than expected for shock breakout through the stellar sur¬ 
face) to be easily detectable. In fact, there have been 
reports of shock breakout detections in Ty pe II SNe on 
the timescale of hours t o several days ([Sc hawinski et al.| 
2008| Ofek et al. 2010 2013 Drout et al. 2014|), which 
require an unusual amount of circumstellar material in 
order to explain the timing of the events, either in the 
form of an extended stel lar envelope inconsisten t with 
stellar evolution models (Schawinski et al. 2008), or a 

(I Ofek 


0-1M© yr 


vigorous, episodic mass-loss ra te of 
et al.||2010 Drout et al.||2014). 

A tew ol the earliest optical discoveries of core-collapse 
SNe have caught the decline of t he “cooling envelope 
phase” following shock breakout jHamuy et al. 1988 


Stritzinger et al. 2002 Dessart et aT 2008 T Observing 
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SN PSl-13arp 

g,r,i Template Image i band detection, 2013 May 9 

host galaxy with SN difference image 



Fig. 1.— PS1 deep gpi,rpi,ipi reference template (left), 1-band image on UT 2013 May 9 (middle), and i-band difference image (right). 
The position of PSl-13arp in the host galaxy is marked with a blue cross. 


the rise to peak of the cooling envelope emission is more 
desirable because it puts direct constraints on the ra¬ 
dius of the progen itor star (Nakar fe SaripOlO Rabinak 


& Waxman 20111. This thermal emission from the ex¬ 
panding shock-heated ejecta lasts only minutes in the 
soft X-rays, but peaks in the UV at hours up to a cou¬ 


ple days after shock breakout (Nakar & Sari 
the time the thermal emission peaks in the op 


2010). By 
icai, other 


radiative processes dominate, i.e. radioactivity (and hy¬ 
drogen recombination in Type II SNe). The timescales 
of the UV peak are more accessible to SN surveys, how¬ 
ever, they still require observing the SN in the UV before 
it brightens in the optical, which can be days later. With 
GALEX we have attempted to circumvent this problem 
by coordinating wide-held monitoring in the UV with 
ground-based monitoring in the optical. 

The large radii of red supergiant (RSG) stars yield 
longer timescales and larger luminosities during shock 
breakout and its immediate aftermath, making them ap¬ 
pealing for observational searches. Stellar evolution cal¬ 
culat ions predict that R SG stars will explode as Type IIP 
SNe (Heger et al.,2003), and Type IIP SNe have been di- 
[tn 


rectly lin ked w ith RSG proge nitors from pre-explosion 
imaging (Smartt et al. 2009). Indeed, the first early 
UV detections of SINe were of two Type IIP SNe, which 
were discovered serendipitously, from overlap between 
the GALEX Deep Imag ing Survey and the CFHT Su- 


et al. 2008 

)• 



The first planned UV and optical joint survey effort, 

the GALEX Time Domain Survey 1 

Gezari et al. 

2013) 

in the NUV, was coordinated wit 

ti the optical 

Pan- 


STARRSl Medium Deep Survey (PS1 MDS), and yielded 
the early detection of the rise to peak of UV emission 
from a Type IIP SN at z = 0.0862 with a cadence of 
2 days. The early UV peak was well fitted by cooling 
envelope emission from the ex plosion o f a red su pergiant 
(RSG) with R = 7OO±2OOR 0 (|Gezari et al.|2010|. While 
the GALEX observations caught the rise to the peak of 
the cooling envelope emission, the shock breakout itself 
was not resolved by the relatively coarse time sampling of 
the observations com pared to the expect ed shock break¬ 
out duration (< lhr). Ganot et al.|(|2014) present a shal¬ 
lower, wide-held joint IV U V GALEX and optical Palo- 
mar Transient Factory (PTF) survey which detected 7 


core-collapse SNe in the UV with a 3 day cadence. Their 
observed detection rate is consistent with core-collapse 
SN rates combined with a fiducial model for cooling en¬ 
velope emission from RSG stars. However, again their 
cadence was not sufficient to catch the short-lived shock 
breakout signal. 

Here we present the results from a high-cadence mon¬ 
itoring program with GALEX which increased the time 
sampling of the GALEX Time Domain Survey and 
GALEX/PTF survey by a factor of two and three, re¬ 
spectively, and has yielded the first potential detection 
of the elusive shock breakout phase in a Type IIP SN, 
albeit with the requirement that the shock breakout oc¬ 
curred outside the progenitor’s surface, in its circumstel- 
lar wind. The paper is organized as follows, in © we 
describe the GALEX and PS1 MDS observations and 
Gemini spectroscopy, in fj3] we compare the UV/optical 
light curve of PSl-13arp to previous early observations 
of Type IIP SNe, in <H] we compare the observed UV 
burst to models for shock breakout from the surface of a 
star, cooling envelope emission, and shock breakout into 
a wind, and in )|5]we conclude with our favored model of 
shock breakout of an RSG with a pre-explosion mass-loss 
rate of « 1O _3 M 0 yr -1 . 

2. OBSERVATIONS 


2.1. Pan-STARRSl Medium Deep Survey 

The Pan-STARRSl system is a high-etendue wide-held 
imaging system, designed for dedicated survey observa¬ 
tions. The system is installed on the peak of Haleakala 
on the island of Maui in the Hawaiian island chain. A 
description of the Pan-STARRS l system , bo th ha rdware 
and software, is provided by Kaiser et al. (|2010). The 
Pan-STARRSl optical design ( fFlodapp et al.||2004 ) uses 
a 1.8 meter diameter //4.4 primary mirror, and a 0.9 m 
secondary, and a 3.3 deg diameter held of v iew. The 
Pan-STARRSl imager (Tonry & Onaka 2009) comprises 
a total of 60 4800 x 4800 pixel detectors, with 10 gm 
pixels that subtend 0.258 arcsec. 

The PS1 observations obtained through a set of hve 
broadband hlters, which we have designated as gpi, 
rpi, ipi, zpi, and yp±. Although the hlter system for 
Pan-STARRSl has much in co mmon with that used in 
previous surveys, such as SDSS ( Abazajian et al.| [2009), 
there are important differences. Further information on 
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the passband shapes is described in Stubbs et al. (2010). 
Photometry is in the “natural” Pan-STARRSl system, 
m = —2.5\og(flux) + m ', with a single zeropoint ad¬ 
justment m' ma de in each band to conform to the AB 
magnitude scale (jTonry et al.||2012). 

Images obtained by the Pan-STARRSl system are 
processed through the Image Processing Pipeline (IPP) 

I Magnier| 2006), on a computer cluster at the Maui High 
Performance Computer Center. The pipeline runs the 
images through a succession of stages, including flat- 
fielding (“de-trending”), a flux-conserving warping to a 
sky-based image plane (with a plate scale of 0.25 arc- 
sec/pixel), masking and artifact removal, and object de¬ 
tection and photometry. Mask and variance arrays are 
carried forward at each stage of the IPP processing. Pho¬ 
tometric and astrometric measure ments pe rformed by 
the IPP system are described in Magnier (2007) and 
Magnier et al. (20081 respectively. 

PS1 has two extragalactic survey modes: the 37 t sur¬ 
vey observes the entire sky above a declination of —30 
deg in 5 bands with a sparse cadence (a few epochs per 
month per filter), and the Medium Deep Survey (MDS) 
is a deeper survey of 10 fields distributed across the sky 
in 5 filters with a daily cadence. This paper uses images 
and photometry from the PS1 MDS. Observations of 3—5 
MD fields are taken each night and the filters are cycled 
through in the pattern: gp\ and rpi in the same night 
(dark time), followed by *pi and Zpi on the subsequent 
second and third night, respectively. Around full moon 
only ypi data are taken. A nightly stack image is cre¬ 
ated using a variance-weighted scheme from 8 dithered 
exposures, allowing for the removal of defects like cosmic 
rays and satellite streaks, and resulting in a 5 a depth of 
~ 23 mag in the <?pi,rpi,zpi, zpi bands, and ~ 22 mag 
in the ypi band in the AB system. 

The nightly MD stacked images are processed through 
a frame-su bt ract ion analysis using the photpipe pipeline 
(Rest et al. 20051, and significant excursions in the dif- 


ference im ages are detected using a modified version of 
DoPHOT (Schechter et al. 1993). Reference templates 
are constructed tor the Ml) fields from two epochs ob¬ 
tained at the beginning of the observing season. Spec¬ 
troscopically confirmed supernovae are later reprocessed 
using deep template images, constructed from all nightly 
stacks with good seeing in the seasons that the SN did 
not explode, with the routine transphot (Rest et al. 20141 
which imp lements the HOT PANTS image differencing 

, then PSF photometry is 


algorithm (Becker et al. 
applied using DAUPHUT ( 57etson||1987 |. The weighted 


average position of the SN is calculated, and forced pho¬ 
tometry is then applied in every difference image at that 
position. The PS1 MDS supernova search has yielded 
hundreds of spectroscopically confir med SNe, includin g 
Type la SNe for cosmology studies (Res t et all[2 014]), 
core-c ollapse SNe, including 76 SNe IIP (Sanders et al. 
20141, a s well as several rare, superluminous SNe at high 


redshift jC homiuk et al.|2011 Berger et al. 
et al.|[20TTTr 


2012 Lunnan 


2.2. GALEX High-Cadence Time Domain Survey 

We conducted an extended survey with GALEX in the 
NUV during December 2012 through April 2013 of 4 PS1 
MD fields, in order to perform a high-cadence (1 day) 
time domain survey, and to complete GALEX coverage 


started wit h the GALEX Time Domain Survey (Gezari 
et al.|2013|) for 9 out of the 10 PS1 MD fields. Th (TEUV 
detector on GALEX was no longer operational after May 
2009. The 1-day cadence of the GALEX observations 
was chosen to better probe the cooling envelope emission 
of RSG stars, which has a characteristic timescale of 2 


days (Nakar & Sari 2010; Rabinak & Waxman 2011), 
and was just barely resolved by the 2-day cadence of 
the GALEX Time Domain Survey for the Type IIP SN 
2010aq (Gezari et al.|2010). The exposure times for each 
GALEX visit averaged 1.1 ks, resulting in a 3cr depth of 
23.5 mag in the NUV in the AB system. Here we present 
the joint GALEX and PS1 discovery of Type IIP SN PS1- 
13arp, for which we catch the onset of a luminous UV 
burst with GALEX that precedes its standard optical 
behavior observed by PS1. 

2.3. Imaging and Spectroscopy of PSl-13arp 

PSl-13arp was identified as an optical transient in PS1 
field MD06 by photpipe on UT 2013 April 26 at R.A. 
12:18:25.108, Dec. +46:37:01.30 (J2000), with a spatial 
offset of 1.577 arcsec East and 0.275 arcsec North of 
its host galaxy (SDSS J121824.98+463700.5). Figure [l] 
shows the PS1 difference image detection in the i bamfo 
as well as a color image of the deep reference template 
with only the host galaxy flux. The SN was also detected 
as a variable source by GALEX during monitoring of the 
field PS_NGC4258_MOS107 at the 8.3cr level on UT 2013 
April 15. The host galaxy is detected in the NUV , and 
measured from the 4 epochs before the SN is detected 
by the GALEX pipeline to have NUV = 21.20 ± 0.05 
mag, measured from a 6 arcsec radius circular aperture, 
and including a correction for the t otal energy enclosed 
of —0.23 mag (Morrissey et al. 2007). Table ?? gives the 
SN photometry from GALEX and PS1. 

We obtained spectroscopy of the UV/optical transient 
with Gemini/GMOS on UT 2013 May 4. We obtained 3 
exposures with a 1.0 arcsec slit and the R400 grating of 
1200 sec each, resulting in a spectral resolution of ~ 4 
A, with a dither along the slit, and a change in grating 
central wavelength, in order to avoid bad columns, and 
to cover the gap in wavelength coverage for the grating. 
The spectrum was overscan-corrected, bias subtracted, 
flat-fielded, wavelength calibrated, and flux calibrated 
using standard IRAF routines. Observations of a stan¬ 
dard star (Feige 34) were used for flux calibration and 
corrections for telluric absorption. The 2D spectrum dis¬ 
plays strong, spatially extended emission lines (Ha, [N 
II] AA6583, 6548, [S II] AA6716, 6731) indicative of star- 
formation in the host galaxy, from which we measure a 
redshift oi z = 0.1665 ± 0.001. 

We fitted the rest-frame host galaxy broad-band spec¬ 
tral energy dis tribution from GALEX FUV and NUV 
(Martin et al. 2005), PS1 gp\, rpi, ipi, zpi, yp\ deep 
stacks with u-band coverage from CFHT (Heinis et al., 
in prep.), and infrare d from WISE at 3.4^m and 4.6/im 
(Wright et al. 2010), with the code CIGALE (Noll et 
al 2009), using version v0.3, implemented in Python. 
CIGALE combines a UV-optical stellar SED with a dust 
component emitting in the infrared, and ensures full en¬ 
ergy balance between the dust absorbed stellar emission 
and its re-emission in the IR. The fit requires both an 
old stellar population with an exponentially declining 
star formation rate with r = 6.7 ± 2 Gyr as well as 
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Fig. 2.— Gemini spectrum of PSl-13arp (smoothed by 2 pixels) at ~ 17 rest-frame days since shock breakout, in comparison to SN 
2006bp at 16 rest-frame days since shock breakout. Arrows show the potential P-Cygni absorption features associated with broad H/3, Fell 
A5169, Nal A5892, and Ha, whose central wavelength are marked with dashed lines. 
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Fig. 3.— Ha profile of PSl-13arp at ~17 rest-frame days since 
shock breakout. The line is broad (FWHM = 7500 zb 1200 km/s, 
Gaussian plotted with green lines) and blueshifted by 2200 zb 100 
km/s. Red dashed line indicates zero velocity defined by the host 
galaxy Ha line. Hashed region indicates the location of the tel¬ 
luric A band feature that has been corrected for in the line profile. 
The profile is mostly symmetric, and does not yet demonstrate the 
strong P-Cygni absorption that develops in Type IIP SNe typically 
3 weeks after explosion. 


a recent stellar population with a moderate star for- 
mation rate of 3.6 ± 1.2M 0 yr _1 , with a total stel¬ 
lar mass of log(M/M 0 ) = 10.4 ± 0.1. Note that this 
mass is consistent with the gas-phase metallicity that 
we estimate from the ratio of the strong nebular emis¬ 
sion lines in the host galaxy spectrum. We measure 
log([NII]\6584/Ha) = —0.37±0.06, wh ich corresponds 
to a metallicity of log {O/H) + 12 ss 9 (Kewley & Do- 


pita 2002), and is in good agreement with the mass- 
met a llicrfy relation for star-forming galaxies ( Tremonti| 
et al.[2004 |. 


We extract a ID SN spectrum by first extracting a ID 
spectrum using a 5-pixel aperture centered on the SN 
(including SN and galaxy light), and then extract a ID 
spectrum with an aperture centered on a region 5 pixels 
away from the SN center, with only galaxy light. We iso¬ 
late the SN spectrum by subtracting the ID galaxy spec¬ 
trum from the ID SN+galaxy spectrum. We coadded the 
first two exposures, which had the highest SN signal. We 
show the coadded SN spectrum in Figure [2j in compar¬ 
ison to the spectrum of SN 2 006bp at a simil ar phase 

i. We de- 
IIA5169, 


since shock breakout from Quimby et al. 
tect the P-Cygni absorption features of 
Nal A5892, and Ha. We fit the strong broad Ha line with 
a Gaussian with a FWHM = 7500 ± 1200 km/s (shown 
in Figure pi). The line appears mostly symmetric, and 
blue-shiftecf by 2200 ±100 km/s from the galaxy’s nar¬ 
row Ha line. We only see weak evidence of blue-shifted 
absorption, as one would expect in a P-Cygni profile. 
The weak blueward absorption is consistent with a young 
Type IIP SN, which typically do not develop a strong P- 
Cygni Balmer line shape until ~ 3 weeks after explosion 
(Quimby et al.||2007 Dessart et al. 2008). However, the 
weak P-Cygni absorption could also reveal unusual con¬ 
ditions in the SN envelope or circumstellar environment. 
Unfortunately, without a later epoch of spectroscopy, we 
can only speculate as to whether PSl-13arp would have 
developed a stronger P-Cygni line profiles at late times. 

3. COMPARISON TO TYPE IIP SNe 

In Figure [4] we plot the GALEX NUV plus PS1 gpi, 
rpi, *pi, zpiiight curve of PSl-13arp. The global trends 
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Fig. 4.— GALEX NJJV and PS1 <?pi, ^pi, ipi, £pi light curve of PSl-13arp with lcr error bars, and arrows indicating 3cr upper limits. 
Dotted magenta line traces the ~ 1 day UV burst. Thick black line indicates the time of the Gemini GMOS-N spectrum. Thick solid lines 
show the CMFGEN model for SN 2006bp (modified to the rest frame of SNLS-04D2dc at z = 0.1854, close in redshift to PSl-13arp) in the 
NUV, g , r, i, and z bands from Dessart et al. (2008) scaled by a factor of (1 + £i 3 arp) in the time axis to account for time dilation. We use 
a minimized y 2 fit to the 4 optical bands in order to determine the time of shock breakout (t s bo) °f MJD 56396.49 =b 0.40 in PSl-13arp, 
or 10 d=8 rest-frame hours after the last NUV non-detection, and a vertical offset that is applied to all bands, of —0.058 =b 0.036 mag. 


in the light curve are a ~ 1-day NUV burst, a rapid gp\- 
band decline of ~ 2.5 mag/(100 day), and an over 50 
day plateau in the gpiand ipibands. The optical light 
curve is characteristic of a normal Type IIP SN, and is 
well fitted wi th the spectral synth esis model for Type IIP 


SN 2006bp ( |Dessart et al.| 2008[ ) K-corrected to a simi¬ 
lar redshift to P Sl-l.'iarp (t h e model for SNLS-04D2dc at 
z = 0.1854 from|Gezari et al.|(|2008|)) (shown in Figure[i|. 
The excellent agreement of the colors in the plateau im¬ 
ply a very small amount of internal extinction [Ay <0.1 
mag). We applied a scaling of (1 + 2 13ar p) in the time 
axis to account for time dilation, and we used a mini¬ 
mized x 2 fit to the 4 optical bands determine a vertical 
offset factor —0.058 ± 0.036 applied to all bands, as well 
as a horizontal shift, inferring a time of shock breakout 
for PSl-13arp of MJD 56396.49 ± 0.40. This explosion 
time, determined entirely from the optical bands, is con¬ 
sistent with the timing of the NUV burst, and occurs 
10 ± 8 rest-frame hours after the last NUV upper limit. 
It is also in excellent agreement with the time of explo¬ 
sion determined from Baye sian parametri c model fitting 
to the optical light curve in Sanders et al. (20141 of MJD 
56396.6 ± 0.5. The absolute magnitude ot th e r-band 


plate au at 50 days of M p , r = —17.54±0.12 mag (Sanders 
et al.|[2014 ) is 0.2 mag fainter than SN 2006bp, and well 


within the range of normal Type IIP luminosities. 

We calculate the peak absolute magnitude of the NUV 
burst to be, 


M N uv = m ~ DM — ^Gai = —18.0 mag, 

where we use m = 21.6 mag, dp = 800 Mpc, DM = 39.5 
mag, and Galactic extinction Aq^ = 8.2 E[B—V) = 0.13 


mag, where E[B — V) = 0.016 mag. Note that for a hot 
blackbody temperature of « 10 5 K, the GALEX NUV 
effective wavelength, A res t = 2316A/(1 + z) = 1985A is 
on the Rayleigh-Jeans tail, and the K correction could 
be as large as K[z) ~ — 7.51og(l + z) = —0.5 mag. The 
peak UV magnitude is consistent with the theoretical 


(Gezari et al. 

2008 Schawinski 

et al. 

2008 

Nakar & 

Sari 2010 Ra 

binak & Waxrnan 

2011). 

and is 

brighter 

than the post shock breakout cooling peak that has 
definitively been observed in SNe IIP, SNLS-04D2ck 

with NUV = 

— 17.0 mag ( 

Gezari et al. 

2008) and SN 


In Figure [5] we show a comparison ot U V light curves 


the UV, 

SN SNLS-04D2dc 

(Gezari et al. 

2008 

) and SN 

2010aq 

Gezari et al. 

2010) 

The comparison 

SNe have 


time-dilation. While the optical behavior of the SNe is 
quite similar, there is a clear excess again in the UV 
burst of PSl-13arp, implying a different source of energy 
from the cooling envelope model fitted successfully to 
these comparison SNe. 


4. COMPARISON TO MODELS 
4.1. Shock Breakout in the Stellar Envelope 

While a burst of neutrinos marks the time of core¬ 
collapse, the supernova shock breakout is the first elec¬ 
tromagnetic signature of a SN explosion. The super¬ 
nova shock “breaks out” when the radiative diffusion 
timescale, 
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Fig. 5.— Type II SNe with early UV photometry. Comparison with SN 2010aq (left) and SNLS-04D2dc (right). A constant offset of +1.47 
mag has been applied to the SN 2010aq data from Gezari et al. (2010) and —0.27 mag to SNLS-04D2dc data from Gezari et al. (2008), 
plotted with solid thick lines. For clarity, we do not plot the GALEX FUV light curve of SNLS-04D2dc, which is lower signal-to-noise 
than the NUV light curve. 


3 {SR) 2 _ 3k P (SR) 2 _ 3 tH p 
A c c c 


2008 

Schawinski et al. 

2008 Gezari et al. 2010 Ganot 

et al. 

20141). In the following sections, we address both 


scenarios. 


where SR is the depth of the shock, A = 1/np is the pho¬ 
ton’s mean free path, k is the opacity, r is the optical 
depth, and H p is the atmospheric scale height, is shorter 


than the shock travel time, t s = 5R/v s (Arnett 1996). 
For a standard RSG envelope density and radial struc¬ 
ture, the duration of the radiative precursor is expected 
to be ~ 1000 sec (Klein & Chevalier 1978 Matzner & 


McKee|1999 ). However, the shock breakout pulse will be 


prolonged up to ~ 1 hr due to light-trave l time effects in 


in the presence of a dense circumstellar wind (Ofck et al. 

2010 

Chevalier & Irwin||2011 

Moriya et al. 

2011). 


PSl-13arp are a factor of greater than ~ 50 longer than 
both the radiative diffusion time at shock breakout from 


Matzner & McKee (1999): 


Uo = 860fc- 0 - 58 (p 1 /p*)-°- 28 K+' 79 M 1 0 5 21 i?| 0 1 0 6 sec, 

where R 500 is the stellar radius in units of 500 Rq, M w 
is the stellar mass in units of 10 Mg, and E$\ is the 
energy of the explosion in units of 10 51 erg s” 1 , and the 
light-crossing time of a RSG: 


ti c = R*/c = U60i?5oo sec. 

Even assuming an envelope density a factor of 1000 
less dense than the standard stellar structure models 


of Wooslcy & Hcger (2007), as was done in Schawinski 
et al. pOOS ), one can increase Gbo only by a factor of 
T. However, such an extended, tenuous structure is 
better attributed to a circumstellar shell of material 


rathe r than the stellar envelope itself (Falk & Arnett 
1977). The long duration of the UV burst suggests that 
either the shock is ’’breaking out” in a wind around the 
progenitor star, or it is associated with the post-shock- 
breakout cooling envelope emission previo usly detected 
in the UV in Type IIP SNe with GALEX (Gezari et al. 


4.2. Cooling Envelope Emission 

In the cooling envelope scenario, the GALEX obser¬ 
vations did not catch the brief radiative precursor from 
shock breakout, but rather the UV burst is from the adi- 
abatically expanding shock-heated SN ejecta. As the hot 
envelope expands and cools, the peak energy of the ther¬ 
mal emission shifts from the soft X-rays into the UV 
when 3 kT = hv^uvi^ + z), where vnuv = 2316 A(or 
T = 2.09 eV), causing a peak in the NUV light curve. 
The rate of cooling depends sensitively on the radius of 
the progenitor, since the smaller the progenitor, the more 
energy is lost to adiabatic expansion. Using the analyt¬ 
ical model for the cooling envelope emission for a RSG, 
the l argest likely Type II SN progenitor, from |Nakar fc| 
Sari (2010), this occurs at a time : 


Vak.Nuv = 1.9M^- 23 R 0 50 6 0 s E° 51 20 d a ys. 

The rise-time in PSl-13arp is constrained from the 
time between the last NUV non-detection and the first 
NUV detection to be < 0.87 days in the SN rest- 
frame, or from the fitted time of shock breakout from 
f|3] to be < 0.5 ± 0.3 d. In order for this UV rise¬ 
time to be compatible with PSl-13arp, it would require 
M r 5 °' 23 ^ 500 8£; 5 i 20 < 0.46. The UV peak can occur much 
sooner for a smaller star, such as a BSG (f pea k ~ 0.54 
days). However, the late-time optical light curve of PS1- 
13arp is incompatible with such a compact progenitor. 
For a BSG, with a typical radius of a factor of 10 smaller 
than a RSG, on e would expect a redder and shorter - 
duration plateau (|Kleiser et al.|2011||Dessart et al.|2013|), 
which is inconsistent with the excellent fit of the plateau 
of SN PS-13arp to the SN 2006bp model. 

Figure[6]shows the evolution of the UV/optical spectral 
energy distribution (SED) of PSl-13arp over time, with 
NUV fluxes linearly interpolated to the time of the PS1 
observations. Unfortunately, we do not have a strong 
constraint on the temperature of the UV burst. We can 
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Fig. 6.— Interpolated spectral energy distributions (SEDs) of 
PSl-13arp with time, and their least-squares fitted blackbody tem¬ 
peratures. The time of each epoch is labeled in rest-frame days 
since the estimated time of shock breakout of MJD 56396.5 from 
Figure [4] The temperature of the earliest epoch is constrained by 
a 2 -bana upper limit to be T > 2.2 x 10 4 K. The SEDs for 1.7 
d, 2.5 d, and 4.1 d, have been divided by factors of 2, 5, and 10 
respectively, for clarity. 


only put a lower limit from the 2 -band non-detection 
of T > 2.2 x 10 4 K. Given that the next epoch, at a 
phase (</>) of 1.7 rest-frame days since our estimated time 
of shock breakout, is hotter, with a least-squares fitted 
temperature of T = (2.8 ± 0.3) x 10 4 K, we adopt this 
as our lower-limit to the temperature of the UV burst at 
4> = 0.8d, since the SN should be cooling with time. The 
NUV flux density, corrected for Galactic extinction, then 
corresponds to a peak bolometric luminosity of L bo i > 
1.7 x 10 43 erg s _1 , which is an order of magnitude larger 
than expected fo r the co oling envelope at the time of the 
NUV peak from Nakar & Sari (2010): 


^RSG.peak = 2.70 x 10 42 M^ 087 R 500 E^ 96 ergs- 1 . 

The observed NUV absolute magnitude of PSl-13arp 
is also brighter than the corresponding peak absolute 
magnitude of this model of Mnuv ~ —16 mag. For a 
BSG, the luminosity at the time of the NUV peak is 
even lower (TesG^peak ~8x 10 41 erg s _1 ). Figurej7]com- 
pares the early UV/optical light curve of PSl-13arp with 
the RSG SN shock b reakout a nd cooling envelope emis¬ 


sion model from Nakar & Sari (20101), scaled to match the 
late-time SN 2006bp model fit. While the observations 
and model are in good agreement close to the cooling 
envelope model peak and later (t > 2d), there is a clear 
prolonged UV excess above the model at earlier times. 
We conclude that the peak of the observed UV burst is 
more luminous than expected for the cooling envelope 
emission following shock breakout, and thus we now in¬ 
vestigate whether the UV burst could be from the shock 
breakout pulse itself, but prolonged in time due to the 
presence of a circumstellar wind. 


4.3. Shock Breakout in a Circumstellar Wind 

The ’’shock breakout” is defined as when the radiation 
from the supernova shock is able to escape in front 
of the shock, i.e. when the radiative diffusion time is 


shorter than the shock travel time (see [4.11. For a 
bare star, this occurs in the outer envelope of the star 
(not actually at its surface). However, when the star 
is enshrouded in a dense wind, the shock breakout is 


delayed, and occurs outside the surface of the star, 
in the circumstellar wind. The theoretical literature 


still refers to this 

phenomenon as the shock breakout 

(Chevalier & Irwin 

|2011| |Balberg & Loeb 

2011 

Ofek 

et al. 

shock 

2013), even though the energy deposited from the 
is diffusing out (escaping) into the wind, instead 


of the stellar envelope. A radiation-dominated shock 
in a uniform medium has a characteristic thickness 
in units of optical depth of t s ss c/v s . The shock 
will breakout in the wind if t w > r s , where t w is the 
thickness of the shock traveling through the wind in 
units of optical depth. Here we assume the wind is from 
steady mass loss from the progenitor star M, resulting 
in a density distribution p w = M/4irr 2 v w = Dr~ 2 , 
where D = 5.0 x 10 16 (M o .oi/v w ,io) g cm -1 , and 
Mo.oi = M/(0.01 Mq yr _1 ) and we assume the slow 
wind velocity characteristic of a RSG, = v w /(10 km 


s x ). Chevalier & Irwin (2011) approximate the op¬ 


tical depth tor a shock wave traveling through the 
wind as t w ss R^ 1 nD 1 where Rd is the depth of the 
shock. Hence shock breakout occurs when t w « c/v s , 
such that Rd = nDv s /c. The luminosity of the 
breakout radiation is given by Ef, 0 « E rad Abo, where 
Urad is the radiative energy in the breakout shell, 
and t bo = Rd/v s = kD/c = 6.6fc(M 0 .oiAv.io ) 
days, where t b o is the time of shock b reakout and 


k = k/( 0.34 cm 2 g 4 ). This is given in Chevalier & 
Irwin (2011) as L bo ss -EradAbo, and thus: 


L ho = 6.0xl0 43 fc- 0 - 6 ^i 2 Mf 6 0 - 6 (Mo.oi/«u l ,io) 
with a corresponding temperature of: 


0 2 erg s 1 , 


T rad = 1.1 x 10 5 fc-°- 25 (t d ,i) 


-0.25 


K. 


where t d .i = td/0- day). Given the observed rise-time 
of t r ise A 0.87 d in the SN rest-frame, and the weak de¬ 
pendence of T on td, we assume a shock breakout ra¬ 
diation temperature of T « 10 5 K, which is allowed by 
our observed SED at the UV peak based on the 2 -band 
upper limit during the UV burst. Given our observed 
Galactic-extinction corrected luminosity density in the 
NUV of L v ^ 1+z ) Vo = f„ o 4:Trdl/(l + z) = 6.25 x 10 27 erg 
s _1 Hz -1 , and assuming the radiation from an optically 
thick blackbody, for which L v = (2nh/c 2 )i/ 3 /(e h,y / kT — 
l)47ri? 2 b , one gets i? bb = 5.8 x 10 13 cm and a bolo¬ 
metric luminosity of 2.4 x 10 44 erg s _1 (a factor of 14 
greater than the observed lower limit). The inferred ra¬ 
dius is close to the upp er range o f radii measured for RSG 
(2 —6 x 10 13 cm; Levesque et al. (20051), and potentially 
situated outside the stellar envelope, m the circumstellar 
wind. Given the rise-time of A 0.87 day, this radius then 
corresponds to a shock velocity of v s A 7 x 10 3 km s _1 , 
similar to the velocity width of the detected broad Ha 
line. 


The Chevalier & Irwin (2011) model places two inde¬ 


pendent constraints on the mass-loss rate (M). The first 
is from the timescale of the event. The observed rise time 
corresponds to M A 1.3 x 10~ 3 (v w po/k)td,i Mm yr 1 . 
The duration of th e shock breakout signal in the Cheva- 
lier & Irwin (20111 wind model is on the order of the rise- 
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Fig. 7.— Light curve of PSl-13arp on a logarithmic time scale, in order to comp are the early UV e mission to models for shock breakout 
and the cooling envelope in more detail. Thick lines show the SN IIP model from|Dessart et al.|(|2008|) fitted to PSl-13arp in Figure [d| thin 
solid lines show the numerical model for shock breakout from the stellar surface and its cooling envelope from |Gezari et al.| (;2 008| (G08) 

for a RSG with R = 805 /it. . M _ 1 2.58 and E = 1.2 X 10 51 erg, dotted lines show the analytical surface shock breakout and cooling 

envelope model from|Nakar & Sari] |2010| (NS10) for R = 500-Rq, M = 15 Mg, and E = lx 10 51 erg. The G08 model has been shifted 
vertically to overlap with the DU8 model fit to PSl-13arp at t > 3 d in the NUV. The NS10 model has been shifted vertically to match 
the NUV data point for PSl-13Rp at t ~ 3 d. Note the clear excess in the observed NUV emission above the cooling envelope emission 
components for both models at t ~ 1 — 2 d. 


time to peak. Thus the rise time in PSl-13arp cannot be 
much shorter than the duration of the UV burst of (t, 
0.82 d in the SN rest frame. Thus we favor a mass-loss 
rate closer to 10 _ 3 M Q yr -1 . The observed luminosity 
corresponds to M ~ 1 x 10 -5 (v w ^ w )k~ 3 E^M^ 3 M@ 
yr -1 . Given the strong sensitivity of the peak lumi¬ 
nosity to It si, this could be easily reconciled with the 
mass-loss rate estimated from the timing of the event 
by adopting an explosion energy of U 51 ~ 1.5. Figure 
[ 8 j plots the parameter space for peak bolometric lumi¬ 
nosity versus rise time to peak in the NUV for shock 
breakout in a wind, and cooling envelope emission for 
BSGs and RSGs with a range of = (0.5 — 1.5) x 10 51 
erg, R = 250 - 8501?© for the RSG and R = 25-85 
for the BSG, and M = 10 — 25 M & . The observed peak 
luminosity of PSl-13arp is well above that expected for 
the cooling envelope phase, and is within the range of 
luminosities and timescales for a burst powered by shock 
breakout in a wind. 

The upper limit on Al from the observed rise time is 
consistent with the lack of excess thermal radiation in 


the optical bands at later times. Moriya et al. (2011) 
model the light curves of Type IIP SJNe as a function 
of M , and find that for M > 10 _3 Mg yr -1 , there is an 
observable excess in the UV and optical bands due to the 
conversion of the kinetic energy of the ejecta into thermal 
energy via interaction with the CSM. The mass-loss rates 
in RSG inferred from spectroscopic stu d ies range fro m 
10^ 7 — 10 -4 M yr -1 ( Mauron fc Josselinp iOll [van Loon 


et al.| [2005). Thus while the inferred mass-loss rate from 
the luminosity and duration of the shock breakout is on 
the high end of the observed distribution, these studies 
are sensitive to the time-averaged mass-loss rate, and 
do not sample short-lived, episodic enha nced mass-lo ss 
that may occur in pre supernova-RSGs (Meynet et al. 
2014|), perhaps as a result of pulsation-driven superwinds 


(Yoon & Canticlloj20ld). Unfortunately, we do not have 


sufficient signai-to-noise in the spectrum of PSl-13arp to 
look for a high-velocity absorption feature or ’’notch”, 


by X-rays 

produced by circumstcllar interaction I 

Chugai 

et ah 2007 

)• 



conservative upper limits on the mass-loss rates of pre- 
SN RSGs to < 1(T 6 - lCr 5 M Q yr " 1 from the lack of 
thermal bremsstrahlung emission or synchrotron radia¬ 
tion fr om the interaction of the SN shock wave with the 
CSM (Dwarkadas 2014 Chevalier et al. 20061. Radio 
and X-ray emission from Type UP SIN 1999em was in¬ 
terpreted as the result of circumstcllar int eraction with 


a stel lar wind of M ~ 2 x 10 6 M© yr 1 (Pooley et al. 


2002|. Mass-loss rates should scale with RSG luminos¬ 


ity, and thus with stellar mass. Sanders et al. (20141 


were not able to put a direct constraint on the progenitor 
mass of PSl-13arp from the comparison of the bolomet¬ 
ric light curve to hydrodynamic models, however they did 
infer a ejected nickel mass from a comparison of the late¬ 
time bolometric light curve of PSl-13arp to SN 1987A 
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Fig. 8 .— Parameter space of peak bolometric luminosity vs. the rise time to peak in the NUV covered by models for shock breakout 
in the stellar envelope for a RSG (red star) and BSG (blue star), shock breakout into a wind (grey bars), and shock cooling envelope 
for a BSG (blue dots) and RSG (red dots) with a range of values for M, R, M, and E. The size of the dots are scaled by the relative 
size of R within the range of values of 250 — 850-Rq for RSGs, and 25 — 85 Rq for BSGs. The ranges of the other SN parameters are 
M = 10 — 25Mq, E = (0.5 — 1.5) x 10 51 erg. The wind mass-loss rate ranges from M = (.01 — 5) X lO - 3 M 0 /yr, with some values labeled 
on the diagram. The peak bolometric luminosity, and UV rise time are plotted for PSl-13arp (black dot), SN 2010aq (black square), and 
SN SNLS-04D2dc (black triangle). The peak bolometric luminosity of PSl-13arp is plotted ranging from a lower limit for the observed 
SED with T > 2.8 x 10 4 K, up to the expected wind shock breakout model temperature of T « 10 5 K. 


of Mn; = 0.34^022-^0. This amount of nickel is on 
the high-mass tail of the distribution for their SNe IIP 
sample, and could be indicative of a more massive pro¬ 
genitor RSG for PSl-13arp. Another clue is the weak P- 
Cygni absorption in the Ha profile of PSl-13arp at t ~ 17 
d. Weak P-Cygni absorption has been found in brighter 
than normal Type IIP SNe, and explaine d by lo wer enve- 
lope masses and higher mass-loss rates (Gutierrez et al. 
20141. 


5. CONCLUSIONS 

We present the joint GALEX UV and PS1 optical dis¬ 
covery of a Type IIP SN with a normal optical light 
curve, including the trademark plateau, preceded by a 
UV burst with the peak luminosity expected from shock 
breakout, but with a longer than expected duration. We 
rule out shock cooling emission from a compact BSG 
progenitor, and favor a shock breakout from a RSG with 
a high, but not extreme, amount of pre-SN mass-loss 
(« 1CU 3 Mq yr -1 ). This mass-loss rate is high enough to 
increase the height at which the SN shock breaks out to 
beyond the stellar surface, but not large enough to effect 
the optical light curve or spectrum measured weeks after 
the explosion. In particular, we do not detect narrow 
emission lines that could be powered by interaction with 
a dense CSM (from M ~ O.OIM^ yr -1 to O.IM^ yr - 1 ' 1 


as seen in Type Iln SNe (Turatto et al.||1993| |Chugai & 
Danziger||1994 |Smith et al.||2008 j 
However, we must caution that the theoretical models 


discussed here have yet to be verified in the normal case 
of shock breakout in the stellar envelope. In the case 
of Type IIP SN SNLS-04D2dc, the early GALEX UV 
photometry could be fitted either with a delayed shock 
breakout in a non-standard, extended, low-densit y enve¬ 
lope of a RSG fo llowed by the cooling envelope (Schaw- 
i nski et aT| |2008j), or a numerical simul atio n for_oniy the 


cooling envelope component (Gezari et al. 2008). The 
cadence (~ 1.5 hr) and low signal-to-noise photometry 
(S/N ~ 2 — 4) in this source were not enough to favor 
one model over the other, let alone constrain their de¬ 
tailed parameters. While a wind effectively smears out 
the shock breakout signal, and makes it possible to re¬ 
solve with a coarser (1 day) cadence, as we report here 
for PSl-13arp, this interpretation would have stronger 
footing if the theory were proven in the more common 
case of shock breakout in the SN envelope. 

Rapid-cadence observations on the timescale of min¬ 
utes, particularly in the UV band, are required to re¬ 
solve the putative rapid rise from shock breakout in the 
SN light curve, and its second UV peak due to the cool¬ 
ing envelope. Future observations of shock breakout of 
Type IIP SNe in the UV could be a sensitive probe of 
the mass-loss history of pre-SN RSG_s, which is difficult to 
constrain directly (Meynet et al. 2014). Furthermore, the 
presence of a wind will greatly enhance the ability o f fu¬ 
tur e UV t ime domain surveys, such as ULTRAS AT (Sa- 


giv et al. 
of the S-N 


2014), to directly detect this diagnostic phase 
explosion. 
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